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^ . Abstract. In this paper we estimate the initial He content in about 30% of the Galactic globular clusters (GGCs) 

' from new star counts we have performed on the recently published HST snapshot database of Colour Magnitude 

, Diagrams (Piotto et al. 2002). More specifically, we use the so-called _R-parameter and estimate the He content 

• from a theoretical calibration based on a recently updated set of stellar evolution models. We performed an 

CO ' accurate statistical analysis in order to assess whether GGCs show a statistically significant spread in their 

^«.^ , initial He abundances, and whether there is a correlation with the cluster metallicity. As in previous works on 

' the subject, we do not find any significant dependence of the He abundance on the cluster metallicity; this 

provides an important constraint for models of Galaxy formation and evolution. Apart from GGCs with the 
. bluest Horizontal Branch morphology, the observed spread in the individual helium abundances is statistically 

CL(' compatible with the individual errors. This means that either there is no intrinsic abundance spread among the 

JL I GGCs, or that this is masked by the errors. In the latter case we have estimated a firm la upper limit of 0.019 

^ ■ to the possible intrinsic spread. In case of the GGCs with the bluest Horizontal Branch morphology we detect a 

significant spread towards higher abundances inconsistent with the individual errors; this can be fully explained by 
, additional effects not accounted for in our theoretical calibrations, which do not affect the abundances estimated 

for the clusters with redder Horizontal Branch morphology. In the hypothesis that the intrinsic dispersion on the 
individual He abundances is zero, taking into account the errors on the individual i?-parameter estimates, as well 
as the uncertainties on the cluster metallicity scale and theoretical calibration, we have determined an initial He 
' ^ • abundance mass fraction Yggc ~ 0.250 ± 0.006. This value is in perfect agreement with current estimates based 

on Cosmic Microwave Background radiation analyses and cosmological nucleosynthesis computations. 
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1. Introduction 

The determination of the initial He abundance of Galactic 
Globular Cluster (GGC) stars plays an important role in 
astrophysics, because of its wide-ranging implications. 

First of all, due to the fact that GGC stars are 
the oldest objects in the Galaxy, their initial He abun- 
dance (Iggc , where Y denotes the mass fraction of He) 

Send offprint requests to: Maurizio Salaris 

* Based on observations with the NASA/ESA Hubble Space 
Telescope, obtained at the Space Telescope Science Institute, 
which is operated by AURA, Inc., under NASA contract NAS5- 
26555, and on observations retrieved with the ESO ST-ECF 
Archive. 



should mirror the primordial He abundance (Yp) pro- 
duced during the big bang nucleosynthesis. Secondly, the 
presence (or absence) of gradients in Yggc func- 
tion of the cluster metallicity provides information about 
Galactic nucleosynthesis and evolution. Thirdly, the ex- 
istence of an abundance spread at constant metallicity 
would strenghten the case for He being the so-called 
'second parameter' (beside metallicity) that determines 
the morphology of the Horizontal Branch (HB) in the 
GGC Colour-Magnitude-Diagrams (CMDs - Sandage & 
Wildey 1967). 

Empirical estimates of Yggc 7 are necessarily indirect, 
since He lines are not detectable in GGC star spectra, 
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apart from the case of hot HB objects, whose atmospheres 
are however affected by gravitational setthng and radia- 
tive levitation, which strongly alter their initial chemi- 
cal stratification (see, e.g., Michaud et al. 1983). Yqqq 
estimates take advantage of results from stellar evolu- 
tion theory, which show that the evolution of low mass 
Population II stars is affected by the initial He con- 
tent. The so-called i?-parameter (Iben 1968; Buzzoni et 
al. 1983; Caputo et al. 1987), defined as the number ra- 
tio of HB to Red Giant Branch (RGB) stars brighter than 
the HB level {R = Nhb/Nrgb), is employed to determine 
Yggc- At a given mctallicity, a higher initial Hc-contcnt 
implies a brighter HB and, in turn, a lower value of Njigb 
{Nhb is only slightly affected), with the consequent in- 
crease of R. Other parameters derived from stellar evo- 
lution can also be employed (see, e.g., the discussions in 
Sandquist 2000; Zoccali et al. 2000), but they are better 
suited to determine relative He abundances than absolute 
ones. 

In Cassisi, Salaris & Irwin (2003, hereafter Paper I) 

we have shown how i?-parameter determinations in two 
samples of GGCs (Sandquist 2000, 43 objects; Zoccali 
et al. 2000, 26 objects), coupled with a theoretical cali- 
bration obtained from updated stellar evolution models, 
provide a constant value of IggCj with no trend with 
respect to [Fe/H]. The weighted average of the individ- 
ual cluster abundances provides yGGC=0.244±0.006 or 
i^GGC=0-243±0.006, the negligible difference being due 
to the choice of the observational sample. This value of 
Yggc is in agreement - within the errors - with re- 
cent determinations of the cosmological baryonic mat- 
ter density (fife) from the Cosmic Microwave Background 
(CMB) power spectrum obtained by the BOOMERANG, 
DASI, MAXIMA and WMAP experiments (e.g. Pryke 
et al. 2002; Odman et al. 2003; Sievers et al. 2003; 
Spergel et al. 2003) which, coupled to BBN calcula- 
tions (Buries, NoUett & Turner 2001), do consistently 
provide yp=0.248±0.001. It is important to recall that 
earlier studies (Sandquist 2000; Zoccali ct al. 2000) 
that employed older generations of stellar models, found 
^GGC^0.20, in complete disagreement with the CMB con- 
straint and spectroscopic determinations from HII re- 
gions (1^=0. 234±0. 002 according to Olive, Steigman & 
Skillman 1997; Fp=0.244±0.002 according to Izotov & 
Thuan 1998). 

As for the abundance spread at a given [Fe/H], the 
results of Paper I were contradictory, in the sense that 
ii-parameter values from the Zoccali et al. (2000) sample 
show a spread that is completely consistent with the mea- 
surement errors, whereas the sample by Sandquist (2000) 
shows an intrinsic additional la spread in the He abun- 
dance of ~0.03. 

In this paper we measured new i?-parameter values for 
a sample of 57 GGCs from the database of HST GGC ob- 
servations by Piotto et al. (2002). We use the HST flight 
photometric system (without reddening corrections) in- 
stead of the Johnson photometric bands. This approach 
avoids potential subtle errors in the estimate of the visual 



magnitudes of the HB. In fact, the transformation to the 
Johnson system requires the knowledge of the cluster red- 
dening, therefore the accuracy of the observed Johnson 
magnitudes depends on the accuracy of the cluster red- 
dening estimate. 

Our determinations of i?-parameter values constitute 
the largest existing homogeneous database, encompassing 
about 30% of the total GGC population. It is therefore 
extraordinarily well suited to address in more detail the 
question of the intrinsic spread in the individual GGC He 
abundances, and the ratio of He enrichment with metallic- 
ity enrichment (AF/AZ). With regard to the latter quan- 
tity, its evaluation is extremely important because of the 
connections with stellar yields and mechanism of forma- 
tion of the GGCs, and the determination of their absolute 
ages. In addition, it is relevant to make a comparison with 
the value AY"/ AZ between rs 6 and 12 (depending on the 
sample adopted, and with the assumption that the Oxygen 
abundance is a good tracer of total metallicity) estimated 
by Olive et al. (1997) for extragalactic HII regions, and 
the value ^ 3 estimated by Izotov & Thuan (2003). We 
will also address the hypothesis put forward by D'Antona 
et al. (2002), which suggests the existence of an He en- 
riched stellar component within individual clusters, due 
to chemical pollution by the ejecta of massive asymptotic 
giant branch stars; in particular, we will assess if this 
phenomenon might seriously affect our determination of 

^GGC- 

In Sect. 2 we briefly describe the cluster sample and 
the theoretical models employed. Section 3 presents our 
results for the i?-parameter measurements and the deter- 
mination of Yqqg- a. summary and conclusions follow in 
Sect. 4. 

2. Cluster sample and theoretical stellar models 

We have employed our large photometric database of 74 
GGCs observed in the HST B (F439W) and V (F555W) 
bands with the WFPC2 (Piotto et al. 2002). The observa- 
tions, pre-processing, photometric reduction and calibra- 
tion of the instrumental magnitudes to the HST flight sys- 
tem, as well as the artificial star experiments performed to 
derive the star count completeness are described in Piotto 
et al. (2002). For each cluster we measured a number of 
stars that ranges from a few thousand to « 47000 (in 
case of NGC 6388). For 57 of our clusters we have been 
able to measure R. The empirical definition of R is the 
same as in Zoccali et al. (2000), i.e., the number of RGB 
stars is computed starting from the level of the observed 
F555W magnitude of the Zero Age HB (ZAHB; the lower 
envelope of the observed HB star distribution). The deter- 
mination of the ZAHB level has been discussed in Riello 
et al. (2003) and Recio-Blanco et al. (2004, in prepara- 
tion). The error on the observed R values is computed 
by combining in quadrature the Poisson error associated 
with the measurement of Nhb and Nj^gb, and the contri- 
bution due to the uncertainty on the ZAHB level (which 
translates into an additional uncertainty on A^hg b , hence 
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on R). Usually, this second contribution is much smaller 
than the Poisson error associated with the number counts 
Nhb and Nhqb- As for the GGC [Fe/H] values, we take 
into account current uncertainties on this quantity by em- 
ploying both the Carretta & Gratton (1997 - CG97) and 
Zinn & West (1984 - ZW84) scales (we assumed an error 
of ±0.15 dex on the individual abundances). The obser- 
vational data relevant to our analysis are summarized in 
Table El 

The theoretical models and isochrones needed to cal- 
ibrate the relationship between R and Y are the same 
as in Paper I, computed for a range of [Fe/H] and Y 
values. A full description of these models will be pre- 
sented in a forthcoming paper (Pietrinferni et al. 2004, 
submitted to ApJ) however, the basic physical ingredi- 
ents relevant for the computation of the theoretical R 
values have already been discussed in Paper I. Here, 
we briefly recall that our models have been computed 
with the a-enhanced metal distribution given in Salaris 
k Weiss (1998), with <[Q;/Fe]>=0.4. The radiative opaci- 
ties (Alexander & Ferguson 1994; Iglesias & Rogers 1996) 
have been computed specifically for our adopted metal 
mixture, and the electron conduction opacities are from 
Potekhin (1999). The nuclear reaction rates from the 
NACRE database (Angulo et al. 1999) have been em- 
ployed, with the exception of the ^^C(q!, 7)"'^^0 reaction. 
For this reaction we employ the more accurate recent de- 
termination by Kunz et al. (2002), based on 7 angular dis- 
tribution measurements of "'^^C(a, 7)^^0 and a consistent 
R-matrix analysis of the process. The claimed relative un- 
certainty of this new rate is half of the uncertainty quoted 
in previous determinations. The relevant energy loss rates 
from plasma-neutrino processes have been taken from Haft 
et al. (1994). We have used the equation of state by Irwin 
et al. (2004, in preparation)^ computed for our actual a- 
enhanced metal mixture. 

Table n] compares, as an example, the adopted ZAHB 
level and corresponding He core mass, HB and RGB evo- 
lutionary times entering our present theoretical calibra- 
tion of the i?-parameter, with the calibration by Zoccali 
et al. (2000), for y =0.245 and two selected metallicities. 
The data for the Zoccali et al. (2000) model have been ob- 
tained by interpolating within their computed grid, since 
they did not specifically compute models for F=0.245. By 
simply comparing the evolutionary timescales that enter 
the R calibration one notices a reduction of ~20% of the 
value of R when passing from Zoccali et al. (2000) mod- 
els to our new ones, as discussed in Paper I. As a further 
comparison to highlight the change in the theoretical cali- 
bration of R due to the improvements in stellar models, we 
also display in the same table analogous data from Cassisi 
et al. (1998) computations, for y=0.245. These models 
were a sort of intermediate step between the models used 
in Zoccali et al. (2000) and in Paper I, and have never been 

^ the equation of state code is made publicly avail- 
able at ftp://astroftp .phys.uvic.ca/pub/irwin/eos/code/| 
eos_demoJortran.tar.gz 
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Table 1. Comparison of the ZAHB levels, He core masses 
(in solar mass units), HB evolutionary times (in Myr) 
and RGB evolutionary times at brighter magnitudes than 
the ZAHB (in Myr) between the models used in Zoccali 
et al. (2000), in Paper I, and the models by Cassisi et 
al. (1998), for two selected metallicities (see text for de- 
tails). 

Zoccali et al. (2000) MviZAHB) WjjTe) Irgb 
Z=0.001 054 0.499 102 63~~ 

Z=0.006 0.72 0.493 112 55 

Cassisi et al. (1998) 

Z=0.001 049 0500 85 58~ 

Z=0.006 068 0.493 94 47 

Paper I 

Z=0.001 058 0.495 88 65~ 

Z=0.006 0.74 0.489 93 56 



used to determine Yqgc- A calibration based on Cassisi 
et al. (1998) models reduces the value of R by ~10% at 
Z=0.001 with respect to Zoccali et al. (2000), whereas at 
Z=0.006 the reduction is by only about 2%. 

The bolometric luminosities and effective tempera- 
tures of the models have been converted to the HST 
flight system using the transformations by Origlia & 
Leitherer (2000), based on the model atmospheres by 
Bessel et al. (1998). In Paper I we have also shown how 
our treatment of the semiconvective regions during the He 
burning phase does agree with independent observational 
constraints coming from the R2 parameter measured in 
a sample of GGCs by Sandquist (2000 - see Paper I for 
details). 

The theoretical values of R are computed by consider- 
ing the RGB evolutionary time of the RGB mass populat- 
ing the relevant isochrone, and the HB lifetime of a star 
populating the middle of the RR Lyrae instability strip 
(log(re//)=3.85). This is strictly adequate only for those 
clusters with an HB populated at the RR Lyrae instabil- 
ity strip and redward (increasing total stellar mass) , since 
the HB evolutionary timescales are basically unchanged 
when moving from the instability strip towards the red 
(see Paper I and references therein). However, stars pop- 
ulating the bluer part of the HB do show different evolu- 
tionary times, which increase for decreasing total stellar 
mass (bluer colour) . At the bluest end of a typical HB the 
increase of the HB evolutionary time with respect to the 
RR Lyrae strip counterpart can amount to about 20 % 
(see, e.g., Zoccali et al. 2000). 

2.1. The effect of He enriched subpopulations in 
individual clusters 

Another effect that might bias our determination of ycGC 
is the He enrichment in some GGC proposed by D'Antona 
et al. (2000). More in detail, to explain the observed CNO 
abundance anomalies and the extended blue HB tails in 
some GGCs, D'Antona et al. (2002) suggested the exis- 
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Fig. 1. Observed i?-parameter values as a function of the 
cluster [Fe/H] for the two adopted metallicity scales. The 
theoretical calibration for F=0.245 and a cluster age of 13 
Gyr is also shown. 



tence of an He enriched stellar component within the clus- 
ters. This enrichment is due to chemical pollution by the 
ejecta of massive asymptotic giant branch stars belong- 
ing to the cluster; this chemical pollution may explain the 
mentioned CNO anomalies and HB colours. In Riello et 
al. (2003) we showed that this supposed enrichment does 
not alter the brightness of the RGB bump in the lumi- 
nosity function of GGCs. Here we estimate the effect on 
the He abundance estimated through the R parameter. 
As in Riello et al. (2003), and similarly to D'Antona et 
al. (2002), we considered a cluster stellar population com- 
posed of 64% of stars with a given initial He abundance as 
adopted in our models - here denoted as 'normal' - , and 
36% of stars with an initial Y randomly distributed be- 
tween the 'normal' value and a mass fraction 0.06 larger. 
In general, for any morphology of the HB, the level of 
the HB that enters the definition of R is still determined 
by the 64% of objects with 'normal' abundance; Nhb for 
the whole population tends to increase, whereas Nugb 
decreases, increasing the helium content estimate. When 
using our present standard calibration of i? as a function of 
Y in He-enhanced clusters, we would overestimate the ini- 
tial He mass fraction by at most 0.005-0.006 for clusters 
where the HB instability strip and the region redwards 
of this strip are populated, whereas for blue clusters the 
effect is an overestimate of at most 0.01-0.02. Most prob- 
ably the degree of He enrichment (if it is real) will vary 
from cluster to cluster (for the clusters affected by this 
phenomenon), and the real overestimate will be between 
zero and the values quoted above. 



3. The value of Iggc 

Figure^displays the observed R values as a function of the 
cluster metallicities, on both the CG97 and ZW84 metal- 
licity scales. The solid line denotes the theoretical values 
computed for an age of 13 Gyr and an initial y=0.245. As 
already discussed in Paper I, theoretical models predict an 
almost constant trend of R with [Fe/H], with a disconti- 
nuity located between [Fe/H] 1.15 and [Fe/H] 0.85. 

The abrupt increase of R in this metallicity range is due to 
the fact that the RGB bump, previously located at bright- 
nesses larger than the ZAHB, moves below the ZAHB level 
with increasing metallicity, thus causing a decrease in the 
number of RGB stars brighter than the ZAHB (see, e.g., 
Salaris et al. 2002). As explained in Paper I, it is only in 
this metallicity range that the theoretical calibration of R 
as a function oiY does depend on the assumed age, due to 
the fact that the RGB bump brightness is also a function 
of age. The net effect is to decrease the theoretical value 
of i? at a given Y when the age decreases. 

Outside this narrow metallicity range a variation of the 
age between, i.e. 14 Gyr and 8 Gyr, produces a change in 
the inferred value of Y by less than 0.001. This is why 
in this case we can safely neglect the precise individual 
values of the cluster ages in this metallicity range, and 
just assume a common reasonable age. Just to show the 
consistency of our GGC isochrone ages with independent 
estimates of the age of the universe (e.g. the WMAP re- 
sults) , we have determined the age of three clusters span- 
ning approximately the entire [Fe/H] range of the GGC 
system, i.e., 47 Tuc, M 3 and M 15; if we consider, e.g., 
the ZW84 metallicity scale, these clusters span the range 
between [Fe/H]=— 2.10 and —0.7. In more detail, we have 
compared the accurate empirical estimate of the pa- 
rameter values (e.g. the V magnitude difference between 
the Turn Off and ZAHB level, which is a strong function 
of the cluster age and weakly affected by metallicity) pro- 
vided by Rosenberg et al. (1999), with the corresponding 
values obtained from our isochrones (transformed to the V 
Johnson photometric system using bolometric corrections 
to the V band homogeneous with the colour transforma- 
tions discussed before). We obtained ages of 11.9±1.2 Gyr, 
11.5±0.6 Gyr, 11.3±1.1 Gyr for, respectively, 47 Tuc, M 3 
and M 15. These ages are consistent with a common value 
between 11 and 12 Gyr, and consistent also with an age 
of the universe of 13.7 Gyr as estimated from the WMAP 
results. Using the CG97 metallicity scale in place of the 
ZW84 one has only a small effect on these ages, of less 
than 1 Gyr. We reiterate again that the exact value of 
the age assumed for the clusters in our analysis does not 
influence at all our results about Yqqc- 

3.1. Estimate ofYcGC using the ZW84 metallicity 
scale 

We have first considered the ZW84 metallicity scale. The 
individual He abundances with associated Icr errors are 
displayed in the lower panel of Fig. El There are 6 clusters 
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Table 2. Cluster data. The columns display, respectively, cluster name, measured value of R with associated error, 
[Fe/H] on the ZW84 and CG97 scale, and value of HBtype (see text for details) 



Cluster 




R 


(7 


{R) [Fe/HJzw84 


[Fe/H 


1 

JCG97 




KJ4499 


1, 


.351 


0. 


.302 


-1, 


.50 


— 1, 


.27 


0, 


,11 


NGC 104 


1 


.607 


0. 


.157 


-0 


.71 


— 


.70 


-0, 


.99 


NGC 362 


1 


.358 


0, 


,193 


-1, 


.27 


— 1 


.15 


-0, 


,87 


NGC 1261 


1 


.208 


0, 


,170 


-1, 


.31 


— 1 


.10 


-0, 


,71 


NGC 1851 


1 


.457 


0, 


,163 


-1, 


.36 


— 1 


.14 


-0, 


,36 


AT /~i -i r\r\ A 

NGC 1904 


2, 


.055 


0, 


,251 


-1, 


.69 


— 1 


.37 


0, 


,89 


NGC 2808 


1 


.598 


0, 


,139 


-1, 


,37 


— 1 


.15 


-0, 


,49 


NGC 3201 


1 


.136 


0. 


,319 


-1, 


,37 


— 1 


.15 


0, 


,08 


NGC 4147 


1 


.767 


0. 


,408 


-1, 


.80 


— 1, 


.59 


0, 


,55 


NGC 4372 


1 


.111 


0. 


,366 


-2, 


.08 


— 1, 


.94 


1, 


,00 


NGC 4590 


0, 


.854 


0. 


,200 


-2, 


.09 


— 1 


.99 


0, 


,17 


NGC 4833 


2, 


.189 


0. 


,405 


-1, 


.86 


— 1 


.58 


0, 


,93 


NGC 5024 


1 


.477 


0. 


,159 


-2, 


.04 


— 1, 


.89 


0, 


,81 


NGC 5634 


1, 


.433 


0. 


,184 


-1, 


.82 


— 1, 


.61 


0, 


,91 


NGC 5694 


1, 


.537 


0. 


,167 


-1, 


.92 


— 1, 


.73 


1, 


,00 


NGC 5824 


1, 


.415 


0. 


,092 


-1, 


.87 


— 1, 


.67 


0, 


,79 


NGC 5904 


1, 


.199 


0. 


,146 


-1, 


.40 


— 1, 


.11 


0, 


,31 


NGC 5927 


1, 


.561 


0. 


,199 


-0, 


.30 


— 0, 


.62 


-1, 


,00 


NGC 5946 


1, 


.321 


0. 


,171 


-1, 


.37 


— 1, 


.15 


0, 


,71 


NGC 5986 


1 


.423 


0. 


,144 


-1, 


.67 


— 1 


.44 


0, 


,97 


NGC 6093 


1 


.031 


0. 


,144 


-1, 


.67 


— 1 


.44 


0, 


,97 


NGC 6139 


1 


.244 


0. 


121 


-1, 


,65 


— 1 


.42 


0, 


,91 


A T y^ 1 y^ 1 -1 -1 

NGC 6171 


1 


.560 


0. 


,432 


-0, 


,99 


-0, 


.87 


-0, 


,73 


NGC 6205 


1 


.719 


0. 


,197 


-1, 


,65 


— 1 


.39 


0, 


,97 


NGC 6218 


1 


.366 


0. 


,292 


-1, 


,61 


— 1 


.37 


0, 


,97 


NGC 6229 


1 


.485 


0. 


,140 


-1, 


,54 


— 1 


.30 


0, 


,24 


NGC 6235 


0, 


.949 


0. 


,218 


-1, 


.40 


— 1 


.17 


0, 


,89 


NGC 6266 


1 


.662 


0. 


,153 


-1, 


.28 


— 1 


.07 


0, 


,32 


NGC 6273 


1 


.554 


0. 


,136 


-1, 


,68 


— 1 


.45 


0, 


,97 


NGC 6284 


1, 


.210 


0. 


,157 


-1, 


.40 


— 1, 


.17 


0, 


,83 


NGC 6287 


1, 


.519 


0. 


,273 


-2, 


.05 


— 1, 


.90 


0, 


,98 


NGC 6293 


1, 


.351 


0. 


,190 


-1, 


.92 


— 1, 


.73 


0, 


,90 


NGC 6304 


1, 


.818 


0. 


,306 


-0, 


.59 


— 0, 


.68 


-1, 


,00 


NGC 6342 


1, 


.771 


0. 


,375 


-0, 


.62 


— 0, 


.69 


-1, 


,00 


NGC 6356 


1, 


.658 


0. 


,153 


-0, 


,62 


-0, 


.69 


-1, 


,00 


NGC 6362 


1 


.429 


0. 


,367 


-1, 


,08 


-0, 


,96 


-0, 


,58 


A T y^l y^i o O O 

NGC 6388 


2, 


.130 


0. 


,121 


-0, 


.74 


-0, 


.74 


-0, 


,67 


NGC 6441 


1, 


.854 


0. 


,113 


-0, 


.59 


-0, 


.68 


-0, 


,78 


NGC 6522 


1 


.183 


0. 


,159 


-1, 


.44 


-1 


.21 


0, 


,71 


NGC 6544 


1 


.500 


0. 


,395 


-1, 


,56 


-1 


,32 


1, 


,00 


NGC 6569 


1 


.547 


0. 


,202 


-0, 


,86 


-0, 


,80 


-0, 


,76 


NGC 6584 


1 


.217 


0. 


,235 


-1, 


,54 


— 1 


,30 


-0, 


,15 
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in the metallicity range affected by the assumed value of 
the cluster age; GGCs in this metallicity range do show 
a large age spread (e.g., Rosenberg et al. 1999; Salaris & 
Weiss 2002), towards values lower than the common age 
of the more metal poor objects. This biases the initial He 
abundance inferred with our assumption of a constant age 
(13 Gyr) towards lower values, unless the individual clus- 
ter age is precisely known and accounted for. These 6 clus- 
ters have been therefore excluded from our analysis; when 
considering the whole remaining sample of 51 clusters, we 
did not find any statistically significant relationship be- 



tween Y and [Fe/H]. We obtain SY/5 [Fe/H]=-0.01±0.01, 
a ratio different from zero by less than its associated 2a 
error. 

Because there is no dependence of Y on [Fe/H], we 
may assume that any spread in the He abundance among 
the clusters is not due to Galactic chemical evolution. To 
investigate the existence of an intrinsic abundance spread 
uncorrelated with [Fe/H] and possible asymmetries (which 
can skew the best estimate of the initial GGC He abun- 
dance when using straight averages or weighted averages) 
of the He abundance distribution around the modal value. 
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Fig. 2. He abundance values as a function of [Fe/H] for 
the 57 clusters analyzed. The two panels refer to the two 
metallicity scales adopted in our analysis. The He abun- 
dances determined for the clusters with [Fe/H] in the 
range comprised between the two vertical dotted lines are 
affected by the choice of the cluster age. The horizon- 
tal dashed line denotes the most probable value for Ygqc 
(see text for details). Filled circles denote clusters with 
HBtype >0.8. 

we have performed the following analysis. Instead of using 
histograms with the cluster He abundance distribution, 
we determined a continuous probability distribution func- 
tion (PDF), following a procedure employed by Sakai et 
al. (1996) for the luminosity function of RGB stars. More 
in detail, the He abundance probability density '^{Y) for a 
generic value of Y is determined by replacing the discretely 
distributed individual abundances by their corresponding 
Gaussian functions (in the hypothesis of Gaussian errors 
on the individual abundances), according to the expres- 
sion 

where Yi and ai are the He abundance and associated er- 
ror of the «th cluster. '^{Y) corresponds to a sum of nor- 
malized Gaussian functions, such that a star with smaller 
error is represented by a Gaussian that is more peaked. If 
Y is the same for all clusters, the PDF should correspond 
to a Gaussian function with dispersion a due only to the 
individual measurement errors. 

Figure 121 displays the PDF for the whole sample, with 
an arbitrary normalization. The peak (i.e. the most prob- 
able value of the He abundance for the entire sample) cor- 
responds to F=0.251, and the la dispersion of the PDF is 



equal to 0.038; these same values are found when consider- 
ing the distribution of the individual cluster abundances. 

To check if this observed dispersion is compatible with 
the cluster individual errors we have performed, as in 
Paper I, the following test. For each individual cluster 
we have calculated a set of synthetic He abundances by 
randomly generating - using a Monte Carlo procedure - 
10000 abundance values, according to a Gaussian distri- 
bution with mean value equal to the observed most prob- 
able value, and a equal to the individual He abundance 
errors. This is repeated for all clusters in the selected sam- 
ple and the 10000 values for each individual clusters are 
combined to produce an "expected" distribution for the 
entire cluster sample, on the assumption that the detected 
He abundance dispersion is not intrinsic, but due just to 
the individual errors, assumed Gaussian. The F-test was 
then applied to determine if this "expected" distribution, 
which has an approximately infinite number of elements, 
displays a variance that is statistically consistent with the 
observed distribution of 51 objects. We state that a Yggc 
range does exist if the probability that the two distribu- 
tions have different variance is larger than 95%. 

When applying the F-test to our data, we find that 
this probability is above 99%, therefore we can formally 
conclude that the observed dispersion is incompatible with 
the individual error bars. Also the shape of the distribu- 
tion is different from Gaussian, as shown by the value of 
the kurtosis, which is equal to 2.8. 

To investigate this matter further, we made a clus- 
ter selection on the basis of their HB colour. We con- 
sidered the value of the ratio HBtype={B-V) /(B+Y+R) 
(Lee et al. 1994), where B, V and R denote the num- 
ber of HB stars, respectively, bluer than the RR Lyrae 
instability strip, inside the strip and redder than the in- 
stability strip. The values of HBtype are from the cata- 
logue by Harris (1996), with the exception of NGC 5634, 
NGC 5946, NGC 6273, NGC 6284, NGC 6388, NGC 6441, 
NGC 6569, NGC 6642, which are not present in Harris' 
catalogue. For these clusters we considered the instability 
strip boundaries of Bono et al. (1995; 1997), and deter- 
mined the number of stars at the blue and red side of the 
strip (appropriately reddened according to the individual 
reddenings given in Harris' catalogue). We then computed 
the value of HBtype by considering V equal to zero, since 
we do not have an estimate of the number of RR Lyrae 
stars for these clusters. In this way we have only an ap- 
proximate upper limit for the value of HBtype, which, as 
we will see in the following, is enough for our purposes. 
By considering various sample selections according to the 
value of HBtype, we have determined that for the 30 ob- 
jects with HBtype <0.8 the F-test gives a probability lower 
than 95% that there is an intrinsic spread in the initial He 
abundances. The observed la spread - equal to 0.027 - is 
therefore compatible with the individual errors. Clusters 
with HBtype in this range (see Fig. |2) still cover the entire 
[Fe/H] interval spanned by the whole sample; again, they 
do not show any statistically significant trend of He abun- 
dance with [Fe/H] {SY/S [Fe/H] =0.002±0.010). The He 
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abundance of this sample is therefore formally constant 
and can be considered to be an estimate of Iggc- The 
PDF for this sample - displayed in Fig. O - provides a 
most probable value y=0.250. It is important to notice 
that the PDF is very symmetric, which is confirmed by 
the fact that the weighted average of the individual abun- 
dances is equal to 0.251±0.003, fully consistent with the 
most probable value provided by the PDF. The kurtosis 
is equal to only —0.09, again in agreement with the fact 
that the PDF can be well approximated by a Gaussian 
function. If we include in this sample with HBtype <0.8 
also NGC 5634, NGC 6273 and NGC 6284, for which the 
computed values of HBtype=0-91, 0.97 and 0.83 are only 
upper limits, the previous results are unchanged. 

This result means that the spread found for the whole 
sample, shown to be inconsistent with the individual mea- 
surement errors, must be due to the 21 clusters with 
HBtype ^0.8; in fact, they show a Icr abundance spread 
of 0.05 around the most probable value of the associated 
PDF F=0.255 (see Fig.EJ. This abundance spread is much 
larger than for clusters with redder HB colours, in spite 
of the fact that there is no trend in the error on the indi- 
vidual He abundances with respect to HBtype. The F-test 
statistics confirms that the abundance spread is inconsis- 
tent with the individual measurement errors at the level 
of more than 99.999%. 

If we do not include NGC 5634, NGC 6273 and 
NGC 6284 in the cluster sample with HBtype >0.8, the in- 
trinsic spread is confirmed with the same confidence level. 

For the 21 clusters with HBtype >0.8 we do not 
find a correlation between their He abundance and 
HBtype', surprisingly, we obtain a correlation with [Fe/H], 
with a slope significant at more than 2a level, i.e., 
(5r/(5[Fe/H]=-0.10±0.04. This correlation is unexpected 
because it does not exist for clusters with redder HB types. 
We will see that it is also absent from the subsample with 
blue HBs if we use the CG97 [Fe/H] scale, and we believe 
it is just an artifact due to the uncertainty on the GGC 
metallicity scale. 

It is therefore confirmed that the large abundance 
spread incompatible with the measurement errors found 
for the whole sample is due to these clusters with the 
bluest HB colours. They are the most affected by the 
possible presence of stellar subpopulations with enhanced 
He, as well as by the dependence of the HB evolution- 
ary timescale on the stellar mass. Both effects tend to in- 
crease the He abundance estimated with our calibration. 
For these blue HB clusters we find that the peak of the 
asymmetric PDF is located at an abundance larger than 
the value obtained for clusters with redder HB, which dis- 
play an abundance spread compatible with the measure- 
ment errors. Figure |31 also shows how the high abundance 
tail of the PDF has more weight than for clusters with red- 
der HB (the value of the kurtosis is equal to 1.0 for this 
subsample of blue HB GGCs). We therefore interpret the 
He abundance spread being incompatible with the mea- 
surement errors found in this group of blue HB clusters 
as the effect of the increasing evolutionary timescale for 
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Fig. 3. Probability Distribution Function (PDF) for the 
He abundance in our GGC sample when using the ZW84 
metallicity scale. The dotted line represents the PDF for 
the 51 clusters with [Fe/H] in the range not affected by the 
cluster age (see text); the solid line represents the subsam- 
ple with HBtype <0.8, while the dashed line represents 
the subsample with HBtype >0.8. The individual PDFs 
are not normalized, but for each value of Y the sum of 
the PDFs for the two subsamples selected on the basis of 
their HB colour is equal to the PDF for the whole sample. 
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Fig. 4. As in Fig. El but for the CG97 [Fe/H] scale (see 
text). 
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very blue HB objects, eventually coupled to the presence 
of He enriched stellar subpopulations. This conclusion is 
reinforced by the results of the next section. 

To avoid these problems with blue HB clusters we will 
therefore consider - as said before - the sample of 21 clus- 
ters with HBtype <0.8 for the best estimate of Yggc- 
Since the observed spread is compatible with the mea- 
surement errors, one possibility - which is what we have 
implicitly assumed in Paper I - is to consider that the 
intrinsic dispersion of the individual Y is zero; this as- 
sumption is fully compatible with the F-test analysis. In 
this case the best estimate of Yggc is the weighted av- 
erage (which also agrees with the mode of the associated 
PDF, as shown before) that provides rGGC=0.251±0.003. 

Another possibility is that there exists an intrinsic 
abundance spread not detected by the F-test, because it is 
masked by the individual errors. We can investigate this 
issue more deeply by assessing what intrinsic dispersion 
in the individual He abundances is compatible with the 
observed distribution of Y . Applying the F-test to the ac- 
tual sample of Y abundances and to synthetic samples 
computed accounting for the actual individual errors plus 
increasingly large abundance spreads (in case of Gaussian 
and uniform intrinsic spread, although of course the math- 
ematical form of the hypothetical abundance spread is un- 
known), we obtained the somewhat expected result that 
an intrinsic la dispersion up to 0.025 is still allowed by 
present data. This means that 0.025 is an upper limit to 
the possible intrinsic spread within our 21 object sample. 

As a further test we selected a subsample of 12 clus- 
ters with individual errors lower than 0.025. These clusters 
span the whole metallicity range covered by the full sam- 
ple, and their individual He abundances do not show any 
trend with metallicity. The corresponding abundance dis- 
tribution shows the same mode and weighted mean as the 
21 object sample, and a reduced la dispersion of 0.021, 
again fully compatible with the individual errors (as as- 
sessed through an F-test). Following the same procedure 
outlined above, we found that the upper limit to the in- 
trinsic la spread allowed by this sample is equal to 0.019. 
It should be realized that this subsample with smaller in- 
dividual errors, still containing a sizable number of objects 
spanning the entire metallicity range, has the same mode 
and weighted mean and shows an He abundance distri- 
bution with a lower dispersion fully compatible with the 
individual errors. This allows one to safely reduce the up- 
per limit of a possible intrinsic la spread to 0.019. 

From our data we can therefore conclude that, if one 
allows for an intrinsic abundance spread, the most proba- 
ble value Iggc is well represented by the weighted aver- 
age given above, since it is practically coincident with the 
mode of the sample PDF. An upper limit of 0.019 can then 
be associated to the la dispersion around Yggc allowed 
by present data. 



3.2. Estimate of Yggc using the CG97 metallicity scale 

We have repeated our abundance analysis considering the 
CG97 metallicity scale (see Fig. In this case, there are 
15 clusters in the metallicity range affected by the cluster 
age, about 25% of the whole sample. In the following we 
will consider only the 42 clusters with [Fe/H]> —0.85 and 
[Fe/H]< —1.15. For these clusters there is again no sta- 
tistically significant trend of He abundance with [Fe/H] 
{SY/S [Fe/H]=-0.017± 0.014). 

The PDF for the He abundance of this sample is dis- 
played in Fig. ^ The most probable value for the He 
abundance is y=0.250 and the la dispersion is equal 
to 0.04. When applying the F-test to this sample, we 
obtain a statistically significant He abundance spread 
at more than 99% confidence level. The shape of the 
distribution is also very different from a Gaussian one, 
with a value for the kurtosis equal to 3.1. As for the 
ZW84 metallicity scale, we find that the 22 clusters with 
HBtype <0.8 do not show a statistically significant He 
abundance spread, nor a trend of He abundance with 
[Fe/H] {SY/S [Fe/H]=-0.002±0.014); for these 22 clusters 
we obtain a most probable value y=0.248. The PDF is 
again highly symmetric and the value of the weighted av- 
erage of the individual abundances, F=0.246±0.004, fully 
agrees with the location of the peak of the PDF; the kur- 
tosis is negligible, equal to —0.03. 

For the 20 clusters with HBtype >0.8, the abundance 
spread is statistically significant (at more than 99.999% 
confidence level) , the most probable value for the He abun- 
dance being equal to y=0.251, and there is no trend of the 
individual abundances with either [Fe/H] or HBtype- As 
for the ZW84 metallicities, the bluest HB clusters have He 
abundances shifted to higher values, and again the high He 
abundance tail of the PDF has substantially more weight 
than for the case of the clusters with redder HB (the value 
of the kurtosis is equal to 1.5). The lack of correlation of 
the individual abundances with [Fe/H] for this subsample 
strongly suggests that the correlation found with ZW84 
metallicities is spurious, due to the actual uncertainty on 
the GGC [Fe/H] scale. 

As for the ZW84 metallicities, we can consider the sub- 
sample of 22 clusters with HBtype < 0.8 to provide the best 
estimate of Yggc- As before, the F-test analysis indicates 
that either there is no intrinsic spread in the individual 
He abundances, or some intrinsic spread is masked by the 
individual measurement errors. In the first case one can 
use the weighted average yGGC=0.246±0.004 as the best 
estimate of Yggc- In the second case, following exactly 
the same steps as for the ZW84 metallicities, wc find that 
the weighted average (again consistent with the mode of 
the PDF) also provides the most probable value for Yggc, 
with a 0.019 upper limit to the la dispersion allowed by 
the present data. 
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4. Summary and conclusions 

In this paper we have presented a detailed study of the ini- 
tial He abundance of the GGC system. We have used the 
most extensive and homogeneous database of i?-parameter 
values to date, and recent updated stellar evolution mod- 
els. As in Paper I, we do not find any statistically signif- 
icant correlation between the individual cluster He abun- 
dances and [Fe/H]. This suggests a very homogeneous 
value of Y for the GGCs, with practically no He abun- 
dance evolution over the entire [Fe/H] spanned by the 
GGC system, at variance with results on the evolution 
of He in extragalactic HII regions, which measure a gra- 
dient of He with Oxygen (hence with [Fe/H]). This result 
on the evolution of He in GGCs strengthens similar con- 
clusions reached in Paper I, and it has to be taken into 
account in studies of Galaxy formation mechanisms. Note 
that Sandquist (2000) also foimd a constant Y for the 
GGC system using two independent indicators of relative 
He abundances, namely the A parameter (difference in 
magnitude between HB and Main Sequence) and the RR 
Lyrae mass-luminosity exponent A. 

As for the existence of an intrinsic spread of the clus- 
ter initial Y, we found that objects with HBtyp^ <0.8 - 
irrespective of their [Fe/H] - do show a remarkably ho- 
mogeneous and symmetric abundance distribution, whose 
spread is fully compatible with the measurement errors. 
One can therefore conclude that either the He abundance 
is the same for all clusters, without any intrinsic disper- 
sion, or that an eventual dispersion is masked by the in- 
dividual errors. In this case we obtain that our data al- 
low a 0.019 upper limit to a possible la intrinsic disper- 
sion. On the other hand, clusters with blue HB colours 
(HBtype >0.8) show a spread in their initial Y incompat- 
ible with the individual error bars, skewed towards val- 
ues higher than those for clusters with redder HB. We 
interpret this result as due to the increase of evolution- 
ary timcscales along the HB phase for progressively bluer 
HB stars, which causes an overestimate of Y when using 
our calibration. This effect may be coupled to the possi- 
ble presence of a subpopulation of He enhanced stars in at 
least some GGCs (see D'Antona ct al. 2002), which also 
leads to an overestimate of Y (at most by 0.01 - 0.02) 
when determined from our calibration. 

As already noticed, for clusters with HBtype <0.8 we 
obtain a distribution of Y values very close to Gaussian 
and a la dispersion consistent with the measurement er- 
rors. This occurrence goes against the possible existence 
of a sizable spread towards higher Y values due to sub- 
populations of He enriched stars, as in blue HB clusters; in 
fact, in case redder HB clusters, the effect of these subpop- 
ulations is at most at the level of ~0.005, and probably 
smaller. It is also reasonable to assume that these sub- 
populations might not be present in all clusters, and the 
net effect on our estimates of Yggc should therefore be 
negligible. An extreme possibility is that all clusters with 
HBtype <0.8 are affected in the same way by this enhance- 



ment, whose maximum effect would be to bias our Iggc 
estimate by ^0.005 towards too high values. 

We conclude by providing a best estimate for Fggc , in- 
cluding the effect of the still uncertain [Fe/H] scale and the 
sources of systematic errors discussed in Paper I (i.e., the 
error on the -'^^C(q!, 7)"'^^0 reaction rate and the method 
adopted to suppress the breathing pulses during the fi- 
nal stages of central He burning). We have considered as 
reference value the Yggc abundance determined with the 
ZW84 [Fe/H] scale. Starting with this reference Iggc we 
have generated a set of 100000 synthetic He-abundance 
values, by applying (through a Monte Carlo simulation) 
to the reference value a set of random and systematic er- 
rors, according to a given probability distribution. In par- 
ticular, random errors have been modeled according to a 
Gaussian distribution with mean value equal to the ref- 
erence one, and la dispersion equal to the correspond- 
ing random error on Yggc- The systematic uncertainties 
due to the choice of the [Fe/H] scale (which causes a de- 
crease of Iggc by 0.005 with respect to the reference 
value), ^^C(a,7)^^0 reaction rate (variation by ±0.008), 
and breathing pulses suppression technique (increase by 
0.003) have been modeled using a uniform distribution 
spanning the appropriate range. 

The mean value for the final synthetic distribution of 
He abundances is yGGc=0.250±0.006. An intrinsic dis- 
persion with a firm la upper limit of 0.019 around this 
value of Yggc is a priori possible given the observational 
errors. This estimate of Yggc is in good agreement with 
the primordial He abundance inferred from the CMB in 
conjunction with Big-Bang nucleosynthesis computations. 
Within the respective la errors this value is also in agree- 
ment with the results from Paper I. 
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